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ABSTRACT 



Aims. We investigate the peculiar properties of a large TeV y-ray flare from Mrk 501 detected during the 2009 multiwavelength 
campaign. 

Methods. We identify the counterpart of the flare in the Fermi LAT telescope data and study its spectral and timing characteristics. 
Results. A strong order-of-magnitude increase of the very-high-energy y-ray flux during the flare was not accompanied by an 
increase in the X-ray flux, so that the flare was one of the "orphan"-type TeV flares observed in BL Lacs. The flare lasted about 
1 month at energies above 10 GeV. The flaring source spectrum in the 10-200 GeV range was very hard, with a photon index 
r = 1.1 ± 0.2, harder than that observed in any other blazar in the y-ray band. No simultaneous flaring activity was detected below 
10 GeV. Different variability properties of the emission below and above 10 GeV indicate the existence of two separate components 
in the spectrum. We investigate possible explanations of the very hard flaring component. We consider, among others, the possibility 
that the flare is produced by an electromagnetic cascade initiated by very-high-energy y-rays in the intergalactic medium. Within such 
an interpretation, peculiar spectral and temporal characteristics of the flare could be explained if the magnetic field in the intergalactic 
medium is of the order of 10"'^ - 10"'* G. 

Key words. Keywords should be given 



1. Introduction 

Variable y-ray emission from BL Lacs is most commonly inter- 
preted as Doppler-boosted inverse Compton (IC) emission by 
high-energy electrons in relativistic jets aligned with the line 
of sight. In the Synchrotron-self-Compton (SSC) model, the 
seed photons for the IC scattering are provided by the lower 
energy (radio-to-X-ray) synchrotron emission produced by the 
same electrons. Since both the synchrotron and IC emissions 
are produced by the same electrons, the variability of the y- 
ray emission is expected to be accompanied by a similar vari- 
ability in the X-ray emission. However, this is not always the 
case. Strong "orphan" y-ray flares in the TeV energy band are 
sometimes observed without acc ompanying X - ray sy nchrotron 
flares (iKrawczvnski et al.L |2004 iDaniel et al.l l2005f) . The na- 
ture of these orphan flares is not clear. Different mechanisms 
have been p roposed to explain these orph an flares in t he context 
of lep tonic (iKusunose & Takaharall2006h or hadronic (iBottcheit 
I2005h models of blazar emission. 

In the framework of the SSC model for the y-ray emission, 
the emission spectrum is expected to be characterized by a pho- 
ton index F > 1 .5 (F is defined as the slope of the differential 
spectrum dNy/dE ~ E^^) in the energy range where electron 
cooling via synchrotron and/or IC energy losses is efficient. The 
limiting value of F = 1 .5 is considered as a lower bound on F in 
the analysis of VHE y-ray flux attenuation thr ough interactions 
with Extraga lactic Background Light (EBL) (lAharonian et al.L 
120061 12007allbl) . At the same time, harder than F = 1.5 spec- 
tra could, in principle, be produced in the VHE y-ray band 
as a result of a bsorption of y-rays in e ^e~ pair production in- 
side the source (lAharonian et al l l2008h . due to the presence of 



low-energy cut-offs or ha r denings in the electr on distribution 
(iKatarzvnski et al.L l2006al: iBottcher et"all l2008h or the forma- 
tion of narrow energy distribu tion of emitting particles in the 
source ( lAharonian et al.L [20021) . The observation of blazars with 
intrinsically hard y-ray spectra, harder than F = 1 .5, would have 
important consequences for EBL studies and for the physical 
models describing the y-ray emission. 

In this letter we consider the spectral and timing character- 
istics of a large y-ray flare of Mrk 501 observed during a multi- 
wavelength campaign in 2009 ( lAbdo et aUiloTlh . The flare was 
remarkable both because it was an "orphan" flare (the factor of 
10 increase in the VHE y-ray flux was not accompanied by a 
similar simultaneous increase in the X-ray flux) and because the 
spectrum of the y-ray emission during the flare was much harder 
than the spectrum of the quiescent emission from the source. We 
show that the spectral and timing characteristics of the source 
during the flare suggest the existence of two distinct components 
in the y-ray emission: a non-flaring one, characterized by a softer 
photon index F ^ 1.8, close to the time-averaged photon index 
of the source, and a flaring one, with an extremely hard spec- 
trum F ^ 1.1. We discuss the possible physical origins of this 
hard flaring component. 



2. Fermi data selection and data analysis 

For our analysis we used the publicly available data collected 
by the Large Area Te lescope (LAT) on board the Fermi satellite 
(lAtwood et al.L |2009|) between August 4, 2008 and March 15, 
201 1. We processed the "PASS 6" data using the Fermi Science 
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Fig.l. Lightcurves of Mrk 501 in the 30-300 GeV (thick red 
data points, errorbars are 68% confidence levels), 3-30 GeV (thin 
cyan data points) and 0.3-3 GeV (dark blue thin data points) 
bands. The vertical dashed line marks the start of th e flare de- 
tected in the TeV band by Veritas dAbdo etal.Ll201lh . 



Toolfl We filtered the entire data set with gtselect and gtmktime 
tools following the recommendations of the Fermi tearrQ and 
retained only events belonging to the class 4, whi ch are most 
confide ntly identified with y-rays, as described by lAbdo et all 
(1201 Oah . To estimate the flux from the photon counts in 30-400 
GeV band we used the gtexposure tool. The sp ectral analysis 
was d one using the unbinned likelihood method (iMattox et al.L 
Il996l) . taking into account all the sources from the 1-st year 
Fermi catalog dAbdo et al.L 1201 Obi) within 15° from Mrk 501 as 
well as an additional bright source at the position of a quasar 
B 1638+3952, not listed in the catalog. 

Timing analysis. The long-term lightcurves of Mrk 501 
for the whole period of observations are shown in Fig. [1] for 
three different energy bands. For the highest energy band 30- 
400 GeV, the lightcurve was produced using the aperture pho- 
tometry method, using time bins with an equal signal-to-noise 
ratio = 3. To estimate the diffuse background level in this 
energy bin we used a nearby source-free region of radius 3 ° cen- 
tered on Ra= 256.21°, Dec= 34.97°. In the lower energy bands 
the binning was homogeneous in time, with time bins of 30 days, 
produced using the likelihood analysis. 

A flare with a flux increase by a factor of ~ 10 was read- 
ily identified in the lightcurve in the 30-400 GeV e nergy band. 
The fla re period includes the Veritas flare reported bv lAbdo et al.l 
(1201 lb . The moment of onset of the Veritas flare is marked in 
Fig.[T|by the vertical dashed line. The flare could not be identi- 
fied in the lower energy lightcurves shown in Fig. [T] In princi- 
ple, if the flux enhancement during the flare was moderate and 
the flare duration was much shorter than 1 month, the flare signal 
would have been "diluted" in the lightcurves binned using the 30 
day binning. We have verified that reducing the size of the time 
bins down to several days (the time scale on which the source 
could be significantly detected at lower energies) does not reveal 
the flare in the lower energy lightcurves. 

The duration of the flare in the VHE band was no t con- 
strained by the Veritas observations (lAbdo et al 1 120 lib , since 
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Fig. 2. Waiting time between subsequent photons in 100-400 
GeV band. The grey shaded interval marks the time of the flare. 



their observations stopped 3 days after the onset of the flare, 
when the source was still in the flaring state. In spite of the 
low statistics of the VHE signal in LAT, the LAT data could 
be used to estimate the duration of the flare. Fig. |2] shows the 
time delay between subsequent photons in the 100-400 GeV 
band. During the flare the time delay between subsequent pho- 
tons dropped below 10 days for about one month. During the 
rest of the Fermi exposure the £ > 100 GeV photon count rate 
from the source was about 1 photon per 100 days. To estimate 
the duration of the flare we divide all the 100-400 GeV photons 
into "flare signal" and "background". The "background" count 
rate is < 10"^ ph/day. The time interval ~ 55 d following MJD 
54953 has 5 photons, four of which we ascribe to the "flare sig- 
nal" (the probability for a background photon to come in the 
~ 50 d interval is 0.5). The ~ 55 d time scale can be consid- 
ered as an estimate of the flare duration. The uncertainty of this 
estimate is by a factor of ~ 2. Indeed, if the flare duration was 
~ 100 d, instead of ~ 50 d, the probability for all four "flare sig- 
nal" photons to come in the first half of the time interval would 
be ~ 9%, so that the 100 d flare duration is excluded at ^ 90% 
confidence level. The waiting time between subsequent photons 
drops down to ~ 2 d at the peak of the flare. Just before the 
peak no photons were detected during the previous 40 d. This 
constrains the rise time of the flare. Assuming that the flux level 
comparable to the peak flux was sustained over a rise time inter- 
val r,., one could estimate an upper limit on Tr- The probability 
that no photons are detected during a time interval ~ 8 d before 
the peak is ~ 2"^* ^ 6%, so that T,. > 10 d is ruled out at > 90% 
confidence level. 

Spectral analysis. The timing behaviour of the source below 
and above 10 GeV is different: the flare appears most promi- 
nently at the highest energies (> 30 GeV) and is not detectable 
at low energies. Taking this into account we have analyzed 
the low- and high-energy part of the LAT energy band sepa- 
rately. We have extracted the spectra in the 0.3-10 GeV and 
10-400 GeV bands for the time period of the flare, between 
Tq - MJD 54953 and Tq + 30 d, corresponding to the highest 
count rate of £ > 100 GeV photons, see Fig.|2] We assume that 
in each energy band the spectrum is well described by a simple 
powerlaw dNy/dE = A (E/EoT^. 
i7] The results of the likelihood analysis of the spectrum in the 
two energy bands are shown in Fig.|3] One can see that the slopes 
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Fig. 3. Fermi spectrum of Mrk 501 for the 30 day period dur- 
ing the flare (shared area, 68% confidence level of normalization 
and slope of the spectr um). Data points i n the TeV band show 
Veritas measurements (lAbdo et al.L l20T l|) of the spectrum dur- 
ing the first 3 days of the flare (solid lines) and in the quiescent 
state (dashed lines). Model curves show emission from electro- 
magnetic cascade initiated by 100 TeV y-rays in the intergalactic 
medium. The solid thick red curve is for intergalactic magnetic 
field 3 X 10 G with a correlation length of 1 Mpc. Dashed 
curves are for B - 10 G (upper curve) and B - 10 G 
(lower curve). Grey data points show the spectrum of the histor- 
ically brightest source s tate in 1997 measured by CAT telescope 
jAharonian et al.L 1 19961) . 



of the spectrum below and above 10 GeV are significantly dif- 
ferent. Below 10 GeV the likelihood analysis gives the slope 
Fie - 1.8+0.2, while above 10 GeV the spectrumis much harder, 
with The = 1.1 +0.2. We note that the highest energy photon from 
the source detected by LAT during the flare period had an energy 
Emax - 188 GeV. Taking this into account we show the power- 
law model for the source spectrum only up to 200 GeV in Fig. 

E 

The fact that most of the highest energy photons from the 
source come in a narrow ti me window around the time of the 
flare was already noticed by lAbdo et aP (1201 ll) . This peculiar- 
ity of the highest energy photon timing was found to lead to the 
hardest source spectrum during the 30 d flare period, with an av- 
erage slope of F = 1.64 + 0.09 in the 0.2-400 GeV band. Our 
analysis of the spectral variability details during the flare sug- 
gests that the flaring activity appears only at the highest energies. 
Taking this into account, we fit the spectrum extracted from nar- 
row energy bins (black thick data points in Fig.|3) with a broken 
powerlaw, instead of the simple powerlaw model. We find that 
the use of the broken powerlaw model improves the quality of 
the fit, with an F-test probability of the chance fit improvement at 
the level of 8%. The best fit value of the break energy is found fit- 
ting the data by the broken powerlaw model is fibreak = 14 GeV. 
The quality of the data is not sufficient for limiting the break en- 
ergy from above, only a lower bound /ibreak > 3 GeV could be 
found (Fig.|4]i. If the break energy is left free during the fitting, 
the range of allowed values of Fhe is wider, so that the lower 
bound on Fhe < 1-5. 

The difference in the spectral slopes and in the timing be- 
haviour at low and high energies point to the presence of two 
separate components in the emission from the source. The "soft" 
component which appears below ^10 GeV is steady while the 
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Fig. 4. Best fit values (cross) and the 68% confidence contour 
for the Ehreak, F parameters of the broken powerlaw model of 
the flare spectrum. 

"hard" component, visible mostly above 30 GeV, is variable dur- 
ing the flare. Two additional arguments which support the hy- 
pothesis for the existence of a slow-variable soft component are 
the closeness of the slope of the spect rum in the 0.3-10 GeV 
band to the time-average slope F ^ 1.85 dAbdo et al.L 1201 Obi) and 
the consistency of a powerlaw extrapolation of the low-energy 
spectrum to the VHE band with the quiescent spectrum of the 
source measured by Veritas and MAGIC (shown by the dashed 
data points in Fig.|3]l. 

3. Discussion: origin of tlie liard y-ray emission 

The photon index of the high-energy component is extremely 
hard. In particular, it is harder than F = 1.5, often adopted as 
a lower bound on the photon index, based on simple SSC mod- 
els for the broad band emission from blazars. Such a limit on the 
photon index of blazars u nderlies the derivation of constr aints on 
the EBL density (see e.g. lAharonian et aT l (l2006ll2007allbl) ).The 
detection of a hard spectral component in the flare of Mrk 501 
has, therefore, important implications for the modeling of blazar 
spectra and the EBL. In the following sections we discuss pos- 
sible mechanism(s) for the production of a very hard spectral 
component. 

3. 1 . Mechanisms intrinsic to the source 

Several models which could explain very hard intrinsic blazar 
spectra in the y-ray band have already been proposed. 

3.1.1. Internal absorption origin 

lAharonian et al.l (l2008h considered the possibility for the for- 
mation of a hard y-ray blazar spectra through the absorption 
of y-rays in the dense soft photon background, with a narrow 
(e.g. thermal) spectrum, y-rays with energies Ey are most effi- 
ciently absorbed through interactions with soft photons of en- 
ergy es - lO-ffi^/lO'^eV]"' eV. Optical depths t ~ 1 - 10 



3 



A. Neronov, D. Semikoz and A.M. Taylor: Very hard gamma-ray emission from a flare of Mrk 501 



would be sufficient to produce a hard spectrum above the en- 
ergy 10[e,/102 eV] GeV. The observation of a hard spectral 
component in Mrk 501 would imply, therefore, the existence 
of a dense field of soft photons with energies e > 100 eV in 
a compact y-ray emission region. If the soft photon field has a 
thermal spectrum, its temperature should he T > 10^ K, much 
higher than the typical temperature of accretion disks in Active 
Galactic Nuclei (AGN) accreting in a radiatively-efficient way. 
At the same time, this temperature range might not be unreal- 
istic for radiatively-inefficient accretion flows, like those found 
in the Fanarof f-Riley 1 radio galaxies such as the nearby radio 
galaxy M 87 ( iNeronov & Aharonianl |2007|) . A potential prob- 
lem for the "absorption feature" explanation of the hard spec- 
trum is that if the soft photon field has a narrow energy distri- 
bution, the suppression of the y-ray signal should become small 

at energies much lower than 10 l^fj/lO" eVj GeV. This means 
that the flaring component of the spectrum should reappear at the 
low-energy end of the LAT energy range, which is not the case. 
This difficulty could, however, be overcome if the soft photon 
distribution has a broad spectrum extending to energies much 
higher than e ^ 100 eV. This might also be typical for the radia- 
tive ly inefficient accretion flows where the emission in the 0.1- 
1 keV range is supposed to be produced via inverse Compton 
(IC) by the accretion flow electrons with a quasi-thermal distri- 
bution dNeronov et al..,2008.) . 



3.1 .2. IC from hard spectrum electron origin 



iKatarzvnski et alj ( l2006ah and iBottcher etal] ( l2008h pro- 
posed that a hard spectra could be produced via the IC mech- 
anism if the spectrum of IC emitting electrons is hard enough. In 
particular, Katarzynski et al. (2006a) consider the possibility of a 
low energy cut-off in the electron spectrum. In this case the spec- 
trum of IC emissio n in the SSC model co uld have a photon index 
as low as F = 1 /3. lBottcher et al.l ( |2008|) consider the possibility 
of a hard spectrum for the case when the seed photons for IC 
scattering are provided by the Cosmic Microwave background 
(CMB) and the electron spectrum has a slope dNe/dE ~ E^^ 
with y < 2. Each of these models carries with it a potential prob- 
lem with relation to source variability; 

On the one hand, in the model of IKatarzvnski et aP (l2006ah . 
an immediate effect of cooling is that an E^^ type spectrum 
forms at energies below the low-energy cut-off, giving rise to 
a F = 1.5 photon spectra. Furthermore, if the synchrotron 
cooling of the highest energy flux occurs in the Klein-Nishina 
regime, even softer, F > 1.5, photon spectra are obtained. The 
cooling times of electrons emitting above 10 GeV are fcooi ~ 
5 (0.3 GIB) days. Thus only fields of ~ 0.1 G would allow the 
hard spectrum IC emission to live for long enough to explain 
the persistence of the hard spectrum during the 30 day flare. 
Furthe rmore, the accelera tion time of 10 GeV electrons in such 
fields (lRiegeretal.Ll2007h . face. ~ 10"^ (0.1 clv^-^f days, leaves 
the question as to the origin of the rise time of the flare unan- 
swered. A possible way out of this difficulty is that the flare 
rise time is determined not by the acceleration time scale, but 
by the dynamical time scale of the system (this would imply, 
however, that there is a pre-existing population of accelerated 
electrons). Cooling times shorter than the synchrotron cooling 
time fcooi could also be achieved if an additional cooling mecha- 
nism is present, such as adiabatic cooling (Finke et al., 2008). 

Similarly, t he pro blem of the cooling time for the model of 
iBottcher et al.l (l2008l) limits its applicability to short flares. For 



this model, the cooling time of electrons emitting ai Ey ~ 10 - 

100 GeV energies is very long, fjc -4x10^ [£'-^,/30 GeV] yr, 
leaving this mechanism inapplicable for the description of a 
short flaring episode. 

SSC type models embedded within the stochastic acceler- 
ation framework, h owever, do not suffer so badly from these 
timescale problems (IKatarzvnski et al.[l2006bh . The steady state 
spectra in this scenario, obtained when stochastic acceleration 
is balanced by radiative losses, can indeed be hard enough to 
explain the observed hard flare of Mrk 501. Furthermore, for 
these models, the longer acceleration time for 10 GeV elec- 
trons in 0.1 G magnetic fields, face. ~ 9 (100 km s^'/^Sa)^ days, 
can more naturally explain the rise-time of the flaring episode. 
The Alfven speeds motivated by this scenario, however, require 
rather dense astrophysi cal environments with n ~ 10^ cm"^ 
dO' Sullivan etal.L 120091) . 

3.1.3. Hadronic models 



Hard intrinsic spectr a in the y-ray range can also arise 
in ha dronic models (Rachen&Meszaros', '1 9981: Aharonian 
20001: iMucke & Pr o theore , 2001; Reim er et al 



2004 



Zacharopoulou et al.L l201 ll) . Variability times as short as 



~ 10 s can originate from the synchrotron cooling time of 
protons in strong ~ 100 G magnetic fields, if the observed 
y-ray emission has proton and/or muon synchrotron origin. A 
hard spectrum, within this framework, can form as the result 
of intrinsic ab sorption of the y-ray flux on soft photon fields 
in the source dZacharopoulou et"an 1201 ll) . in a way similar to 
the leptonic models with intrinsic absorption. Alternaviely, the 
spectrum of accelerated protons mig ht be intrinsically hard and 
have a pileup at the highest energies dXh aroma nl 120001) . Finallv. 
if significant contribution to the source flux is produced through 
muon synchrotron emission, the hard spectrum could arise by 
muons with an energy below a certain critical energy decaying 
before they lose their energy through synch rotron emission 
(lRachen&Meszarosill998l;lReimer et alT.l2004 . 



3.1 .4. Vacuum gap acceleration origin 

The problem of the abse nce of cooling of e lectron s en- 
countered in t he m odels of IKatarzvnski et al.l (l2006al) and 
IBottcher et al.l (l2008l) disappears if one assumes that the ob- 
served y-ray emission comes directly from the electron accel- 
eration region. In this case the hard spectrum of electrons is 
maintained by the balance between the energy gain and en- 
ergy loss rates. Such a possibility is considered e.g. in mod- 
els desc ribing particl e acce l eration in black hole magneto- 
spheres (iBeskin et al.L 19921 lEevinso nl I2000t iNeronov et al 
2001 
12007 



Aharonian & Neronovi 120051: Neronov & Aharonian 
INeronov et al. [ I2009h in which an almost monochromatic 



distribution of particles is maintained at the vacuum gap in 
the magnetosphere. The spectrum of synchrotron and/or cur- 
vature emission from electrons and/or protons accelerated in 
the vacuum gap can be very hard. For example, in the regime 
when the acceleration of particles of mass m and charge e in a 
magnetosphere with magnetic field B is balanced by the syn- 
chrotron loss rate dE/dt ~ kcB ~ -e*B^E^ /m'* = -^'synch 
{k < 1 is acceleration efficiency), typical particle energies are 
E ~ rrP'K^I'^ lie^l^B^I'^). The synchrotron radiation from such par- 
ticles is sharply peaked at an energy Ey - eBE^/m^ ~ Km/e^ ^ 

10~A- l^m/m^j GeV independently of the magnetic field strength. 
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The spectrum of emission below this peak energy is hard (see 
e.g. IVincent & LeBohed (l2010h for an example of the spec- 
trum calculation). In principle, this hardness could be close to 
the limiting case for the synchrotron emission spectrum from a 
monochromatic electron distribution, which is characterized by 
a photon index F - 1/3. If the energy loss mechanism limit- 
ing the maximal proton energies is curvature radiation, the max- 
imal photon energy depends on the magnetic field and on the 

particle energy as Ey ^ 30 [(£'/m)/10'°]^ [M/10''cm]"' GeV, 
where M is the black hole mass which determines the size of the 
acceleration region. The y-ray emission from the vacuum ga p 
is highly anisotropic (see e.g. ( iNeronov & Aharoniani l2007l) V 
Most of the hard spectrum luminosity is directed along the AGN 
jet so that the y-ray flux from the vacuum gap can be comparable 
or higher than the y-ray emission from electrons accelerated in 
the jet. Since the energy deposite d in the accelerated particles i s 
released in the form of the y-rays eronov & Aharoniani l2007l) . 
the y-ray production efficiency in the vacuum gap is extremely 
high. 

3.1 .5. IC from electrons within a cold relativistic flow origin 



Alternatively, a narrow particle energy distribution could be 
maintained in a "cold" relativistic wind with very large bulk 
Lorenz factor, similar to the case of pulsar winds, if such winds 
are ejected by the central supermassive black hole. Hard spec- 
trum y-ray emission could then originate through IC emission, 
produced by the wind propagating through an external radiation 
field, as suggested bv iAharonian et al.. (.20021) . 

3.2. Emission from y-ray induced cascades in fhe 
intergalactic medium 

Very hard y-ray emission spectra can also form through the prop- 
agation of multi-TeV y-rays from their source to Earth. The ab- 
sorption of VHE y-rays through their interactions with the EBL 
initiates the development of electromagnetic cascades. If the en- 
ergy of the primary y-rays is high enough, most of the primary 
source power is transferred into an electromagnetic cascade. The 
energy in the electromagnetic cascade is released in the form of 
detectable secondary cascade y-ray emission with energies be- 
low the energy at which the mean free path of y-rays i s com- 
parable to the distance to the source (lAharonian et al.l Il994t 
lNeronov&Semiko"3 l2007h . This results in the possibile forma- 
tion of a hard y-ray emission spectrum, resulting from the de- 
velopment of a cascade in the intergalactic medium, as was sug- 
gested bv iAharo nian et al. (2002). An example calculation of the 
cascade spectrum in the case of Mrk 501 is shown in Fig. |3] In 
this calculation primary y-rays with energy 100 TeV were in- 
jected at the source. The broad band type spectrum in the 
energy range 0.1 -10 GeV is formed during the propagation of the 
y-ray beam towards the Earth. The observed high-energy cut-off 
at Ecut - 5 TeV sits at the energy for which the photon mean free 
path is about the distance to the source. 

The low energy suppression of the cascade spectrum arises 
due to the influence of intergalactic magnetic fields on the ge- 
ometrical structure of the cascade. If the intervening intergalac- 
tic magnetic field (IGMF) strength were negligible, secondary 
e^e" pairs, as well as cascade y-rays would propagate in the 
same direction as the primary y-ray beam. Small deflections 
of the trajectories of electrons and positrons by weak, but non- 
negligible magnetic fields lead to the time delay of the cascade 



signa l (iPlagal Il995t iMurase et al.l l2008t iNeronov & Semikozl 
l2009h . At energies for which the time delay is much longer 
than the flare duration, the cascade emission is suppressed. 
This time delay increase s with the decrease of the e nergy Ey 
of the cascade photons (INeronov & Semikozl l2009l) Tj, 



10 



1.5 



elay 



[fi^/lO^GeV] [s/IO '^g] d. The time delay sup- 
pression of the cascade flux becomes stronger at lower ener- 
gies. This is illustrated in Fig. [3] where the cascade spectrum 
calculated for magnetic field values of 10"'^ G, 3 x 10"'^ G and 
10 G are shown. The low-energy suppression of the spectrum 
is observed below energies of 100 GeV. Monte Carlo calcula- 
tions which take into account the production spectrum of the 
pairs and of the secondary cascade photons show that at ener- 
gies for which the time delay suppression is important, the spec- 
trum of the cascade emission is hard, with a photon index F ^ 1, 
close to the observed photon index of emission above 10 GeV in 
the Mrk 501 flare. Details of the Monte-Carlo code used for the 
calculation can be found in lTavlor et al.l (1201 ll) . 

The range of IGMF strengths required in the cascade model 
is consistent with the lower bounds on the IGMF which have 
been recently derived from the non-observation of the y-ray 
emission from electromag netic cascades emission from several 
hard spectrum TeV blaza rs Neronov& Vovk (|20I lb:lTavlor et all 
(|20I lh : Ffavecchio et al.1 (l2OI0l) : lDermer et al.l (1201 ll) . The exact 
level of this lower bound depends on the nature of suppression 
of the cascade emissio n signal (time delay or extended nature of 
the cascade signal, see iTavlor et aH (|20I ll) for details) and is in 
the range of 10"^^ - 10"'* G if the IGMF correlation length is 
/Ifi » 1 Mpc. For shorter correlation lengths, the bound strength- 
ens to Z? » 10 G. This cascade explanation of the observed 
hard spectrum of Mrk 501 is valid only if the real IGMF strength 
is close to the previously derived lower bounds. 

The cascade model provides an alternative possible explana- 
tion of the extremely hard spectrum in the 10-200 GeV band. 
The interpretation of the data in terms of this model provides a 
measurement of the magnetic field in the intergalactic medium 
along the line of sight toward Mrk 501. The duration of the 
flare was ~ 1 month. If the magnetic field were weaker than 
B ~ 10"'^ G, the characteristic time delay in the energy range 
Ey ~ 100 GeV would be much shorter than the flare duration and 
no suppression of the y-ray flux would be observed (dotted curve 
in Fig. [3]). On the other hand, if the magnetic field were stronger 
than 10""' G, the type spectrum formed by the time delay 
suppression effect would continue to energies much higher than 
100 GeV. This would place heavy demands on the source TeV lu- 
minosity, requiring a TeV flare luminosity of more than an order 
of magnitude higher than that observed in the 100 GeV energy 
band. Although this possibility can not be fully ruled out (no 
systematic monitoring of the source during the flare was done by 
the ground-based y-ray telescopes), such a scenario would need 
a source flux larger than 3 x 10"'" erg/(cm^ s), in excess of the 
highest historical flux detected from the source during the large 
1997 flare (the spectrum of the 19 97 flare measured by HEGRA 
telescope (lAharonian et al 1 [19961) is shown in grey in Fig.|3). 

An alternative possibility for the cascade contributio n to th e 
Mrk 501 spectrum was considered by Takahashi et al. ( !2012l) . 
Assuming that the observed TeV band flux is intrinsic to the 
source, rather than resulting from the development of an elec- 
tromagnetic cascade, the casca de emission is still ex pected to 
appear in the GeV energy band. lTakahashi et al.l (1201 2l) made an 
attempt to constrain the strength of IGMF along the line of sight 
towards Mrk 501 by searching for time delayed GeV band "pair 
echo" emission following the hard flare discussed in this paper 
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We note, however, that even if the peak flare flux in the TeV 
band was at the level measured by VERITAS, the Fermi data 
are of insufficient quality to provide sensible constraints on the 
IGMF. The reason for this is the high level of the steady-state 
flux from the source in the GeV energy band. From Fig. [Tjone 
can see that no variations in the flux level are detected in the 
0.3-30 GeV band, neither during the flare, nor in the period fol- 
lowing the flare. The source flux in the GeV band is constantly 
at the level of EFe - 2 x 10"" erg/cm^s, which is above that 
ex pected for the pa ir echo flux at any moment of time (see Fig. 1 
of IXakahashi etal] (|2012|) ). Thus, the pair echo signal is not ex- 
pected to be detectable, no matter how weak the magnetic field 
is. At the same time, if the peak flux of the flare was at the level 
indicated by the Fermi measurement in the £ > 100 GeV energy 
band (see Fig. [3]), th e stren gth of the pair echo signal considered 
byljakahashi et al.] (1201 2h would be larger by a factor 3-5, so 
that the cascade emission in the GeV band would be above the 
level of the steady state emission observed from the source. 
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4. Conclusions 

The observation of a very hard spectrum of y-ray emission from 
Mrk 501 during an "orphan" y-ray flare challenges theoretical 
models of y-ray emission from blazars. Several mechanisms for 
the formation of the very hard spectra are able to provide expla- 
nations for the observed source behaviour To clarify the nature 
of the hard flare(s) a systematic analysis of the frequency of oc- 
currence and properties of the hard and/or "orphan" y-ray flares 
in the TeV y-ray loud BL Lacs is needed. 
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